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ABSTRACT 

We present a study of the shapes, orientations, and alignments of Galactic dark matter subhalos in 
the "Via Lactea" simulation of a Milky Way-size ACDM host halo. Whereas isolated dark matter halos 
tend to be prolate, subhalos are predominantly triaxial. Overall subhalos are more spherical than the 
host halo, with minor to major and intermediate to major axis ratios of 0.68 and 0.83, respectively. 
Like isolated halos, subhalos tend to be less spherical in their central regions. The principal axis ratios 
are independent of subhalo mass, when the shapes are measured within a physical scale like rv max . 
the radius of the peak of the circular velocity curve. Subhalos tend to be slightly more spherical closer 
to the host halo center. The spatial distribution of the subhalos traces the prolate shape of the host 
halo when they are selected by the largest Vmax they ever had, i.e. before they experienced strong 
tidal mass loss. The subhalos' orientation is not random: the major axis tends to align with the 
direction towards the host halo center. This alignment disappears for halos beyond 3r2oo and is more 
pronounced when the shapes are measured in the outer regions of the subhalos. The radial alignment 
is preserved during a subhalo's orbit and they become elongated during pericenter passage, indicating 
that the alignment is likely caused by the host halo's tidal forces. These tidal interactions with the 
host halo act to make subhalos rounder over time. 

Subject headings: cosmology: theory - dark matter - galaxies: dwarfs - galaxies: formation - galaxies: 
halos methods: numerical 



1. INTRODUCTION 

Cold dark matter (CDM) halos are not smooth. Early 
numerical simulations of the formation of CDM halos 
lacked sufficient resolution to detect much besides the 
gross features of the mass distribution ([Davis et al.lll985l; 
Frenk et all Il985| iQuinn et all [1986 : Efstat hiou et all 



There is now a large body of work concerning the 
shapes of isolate d CDM halos (iDubinski fc Carlberd 
19911 lKat3 HMD: IWarren et all Il992l: iDubinskJ 11991 



19881 : iFrenk et al. Il988t ). Rapid advances in computa 



tional power and in the efficiency of codes have afforded 
a dramatic reduction in the particle masses and gravi- 
tational softening lengths. Multi-mass techniques (e.g. 
iKatz fc WhiTi 119931 : lBertschingerll2001h allow even fur- 
ther resolution increases in particular areas of interest. 
The resulting high resolution maps of the matter distri- 
bution in CDM halos h ave revealed an astonishing abun- 
danc e of substructure (iKlypin et aTl 119991 : iMoore et al.l 
ll999HDiemand et al.ll2007aL hereafter paper I). In "Via 
Lactea", the most recent and highest resolution CDM 
simulation of a Galaxy scale halo to date, the total num- 
ber of identifiable subhalos has reached ~ 10,000, which 
together account for 5.6% of the total host halo mass (pa- 
per I). Recent progress notwithstanding, the currently 
achievable mass resolution is orders of magnitude above 
the true cutoff in the CDM fluctuation power spectrum, 
and given the observed scale-invariant nature of the sub- 
halo mass function (dN/dh\M oc M _1 ), it is likely that 
the total substructure abundance and mass fraction has 
not yet converged. 

CDM halos are not round. Whereas analytical 
work often treats CDM halos as spherically symmet- 
ric mas s distributions, it h a s been known for some 
time dBarnes fc Efstathiou] 119871 ; lEfstathiou et all 
119881: IFrenk et all Il988t) that in general CDM halo 
shapes show significant departures from sphericality. 



Jing et all 119951: iTormenl Il997t [Thomas et al 



Bullock 2002; riJing fc Sutol l2002l: ISpringel et a! 
Bailin fc Steinmeta [2005; 



Hopkins et all 120051: lAllgood et all 120061 : iBett et all 
2007t lHavashi et al.l 120071 : iMaccio et all l2007t) and 



widespread agreement on a number of findings: CDM 
halos tend to be prolate, they are more spherical in their 
outer regions, more massive halos tend to have smaller 
axis ratios, the moment of inertia is aligned with the 
shape velocity anisotropy tensor, the smallest principal 
axis tends to line up with the a ngular momentum ve ctor, 
etc. (for a recent summary, see lAllgood et al1l2006l ). 

With Via Lactca's unprecedented resolution, we are 
now for the first time able to extend this shape analysis 
to a well resolved subhalo population. This is interesting 
from a theoretical point of view, because of the important 
role that tidal interactions play in shaping the properties 
of the subhalo population. The shapes of subhalos are 
likely to be affected by tidal deformations, and a careful 
analysis of the subhalo shapes might allow us to better 
understand the tidal interactions. 

Knowledge of subhalo shapes is important for several 
types of observational studies as well. Weak lensing is 
becoming a ver y valuable tool for probing cosmolog- 
ical parameters dBrown et all 120031: iBacon et al.1 120031: 
Hamana et al.l 120031: Jarvis et al.l l2003t iHoekstra et al.l 
20061 : iMassev et al.ll2007| ) and constrain ing density pro- 
files of galaxy groups a nd clusters ([Brainer el ' 2004; 
iMandelbaum et al1l2006al ). Any alignment between in- 
trinsic galaxy ellipticity and with other galaxies or the 
local mass density will introduce a bias into the lens- 
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ing sig nal (iHirata fc Sell aid l200l iBridle fe Kingl 12001 . 
In fact, iLee et al.1 (|2005l ) found evidence for an intrinsic 
alignment, with subhalo minor axes preferentially per- 
pendicular to the host halo's major axis, in a numeri- 
cal simulation of four cluster scale halos, resolved with 
~ 2 x 10 6 particles each. They also developed an ana- 
lytical model under the assumption that this alignment 
was caused by the host halo's tidal field, and showed 
that it reproduces the numerical findings well. Observa- 
tional evidence for subhalo alignment has be found in the 
Sloan Digital Sky Survey (SPSS), for galaxies in clusters 
(|Agustsson fc Brainer J20061 iMandelbaum et"afll2006bf ) 
as well as in groups (jFaltenbacher et aLll2007f). Based 
on N-body numerical experiments lCiotti fc Duttal (|1994D 
concluded that the tidal field of a spherical cluster could 
be responsible for the observed radial alignment of clus- 
ter galaxy isophotes. A better understanding of typical 
subhalo shapes and their alignment within the host halo 
might allow a statistical correction of the resulting bias. 

Another observational arena dependent on subhalo 
shapes are stellar kincmatical studies of Local Group 
dwarf galaxies, which are being used to constra i n the 
masses of their DM host halqs (IWilkinson et al.l 12004 



Lokas et al.ll2005tlMufioz et al.ll2005t IWalker et al l 



2006 



Gilmore et al.ll2007t IStrigari et al.ll2007l) . In almost all 
cases the analysis is performed assuming spherical sym- 
metry and often also a constant velocity anisotropy. 
More sophisticated analyses of the dwarf galaxy stellar 
motions will benefit from firm theoretical expectations 
of the intrinsic DM subhalo shapes and should result in 
more realistic models of the underlying DM mass distri- 
bution. 

This paper is organized as follows. In Section [2] we 
briefly outline the technique we employed to determine 
the shape of subhalos. We present the shape parame- 
ters of the Via Lactea host halo in Section [H and move 
on in Section |4] to discuss the dependence of the sub- 
halo's shape parameters on the radius at which they are 
measured, on their distance from the center of the host, 
and on the subhalo's mass. In Section [5] we present re- 
sults concerning the spatial distribution of the subhalos 
within the host halo and the alignment of their ellipsoids 
towards the host halo center. In Section[6]we consider the 
rcdshift evolution of the shapes of a sample of strongly 
tidally affected subhalos. Section [7] contains a discussion 
and summary of our results. 

2. SIMULATION AND SHAPE FINDING METHOD 

The main sample of subhalos analyzed in this work 
stems from the z — output of the "Via Lactea" sim- 
ulation (paper I). This simulation follows the dark mat- 
ter substructure of a Milky- Way-scale halo (M 2 oo = 
1.77 x 10 12 Mq) with 234 million particles. With a parti- 
cle mass of ~ 20, 000 M Q , the simulation resolves around 
10,000 subhalos within ^200 = 388 kpc. The global z = 
properties of the host halo and the substructure pop- 
ulation was presented in paper I and the joint tempo- 
ral evolution of host halo and substructure properties, 
with an empha s is on t idal interactions, was discussed in 
iDiemand et al.l (|2007bl hereafter paper II). Here we focus 
on the shapes of the matter distribution in Via Lactea's 
host halo and subhalo population. 

For the determination of (sub) halo shapes we follow 
the iterative technique outlined in lAllgood et al.l l)2006h . 



This method is based on diagonalizing the weighted mo- 
ment of inertia tensor 



where 



v4 



iVn/q) 2 + (Zn/S) 2 



(1) 



(2) 



is the ellipsoidal distance in the eigenvector coordinate 
system between the (sub)halo's center and the n th par- 
ticle, and q = b/a and s = c/a are the intermediatc- 
to-major and minor-to-major axis ratios, respectively 
(a > b > c). Initially Jy is calculated for all particles 
within a spherical window of radius tq. In subsequent 
iterations we fix a = ro and include only particles with 
r n < rQ. Iteration continues until q and s change by less 
than 10~ 3 . The degree of triaxiality of a halo is quantified 
by the triaxiality parameter introduced by iFranx et al.l 
l(l99lh 



T 



(3) 



A halo is said to be oblate for T < 1/3, triaxial for 
1/3 < T < 2/3, and prolate for T > 2/3. 

For comparison with observational data it is often more 
desirable to constrain the shape of the potential than 
the shape of the density distribution. The potential 
shape has the additional advantage that it is less sensitive 
to local density variations (from substructure for exam- 
ple) and is typically smooth and well approximated by 
conce ntric ellipsoids (jSpringel et al.ll2004t lHavashi et al.l 
12001 . Instead of measuring the potential shape directly, 
or by fitting ellipses to the intersections of isopotential 
surfaces with three orthogonal planes (as advocated by 
ISpringel et alll200l . we diagonalizc the unweighted ki- 
netic energy tensor 



Ki 



(4) 



where the Vi^ n are measured in the rcstframe of the sub- 
halo under consideration. K^ is related to the poten- 
tial energy tensor = ^ xid^/dxj through the tensor 
virial theorem 



Id 2 1, 
2 



^L=2K ij +W ij . 



(5) 



The internal structure of the Via Lactea host halo re- 
mains practically unchanged after z = 1.7 (paper II) and 
we expect d 2 Iij/dt 2 = at z = 0. In that case the eigen- 
vectors of should reflect the principal axes of the 
potential ellipsoid. Note that the assumption of a con- 
stant Iij is likely to fail for subhalos that are being tidally 
stripped, and in this case the principal axes of Kij will 
not necessarily reflect the shape of the potential. Instead 
of the iterative procedure applied for 1^ , we simply diag- 
onalize Kij for all particles within the moment of inertia 
ellipsoid at a given radius. In the following we will refer 
to shapes determined by diagonalization of and K$j 
as the mass and potential shapes, respectively 

3. HOST HALO SHAPE 

At the present epoch the Via Lactea host halo is pro- 
late. We have determined the principal components of 
Iij and Kij as a function of major axis radius. The re- 
sulting shape ellipsoids are depicted in Fig. [1] Neither 
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the shape nor the orientation of the ellipsoids vary much 
out to r2oo- As a function of radius the major axis eigen- 
vectors for both mass and potential shape are aligned 
to within 1°. As expected the potential is significantly 
rounder than the mass distribution, and the major axes 
of the two are aligned to within 5° at all radii. In Fig. [5] 
we plot mass and potential shape parameters q and s and 
triaxiality parameter T at z =0, 0.5, and 1 as a function 
of radius. Recall (from paper II, Figs. 3 and 4) that the 
spherically averaged mass distribution of the Via Lactea 
host halo remains remarkably constant in proper coordi- 
nates after the last major merger at z = 1.7 until today. 
We now see, however, that the host halo does undergo 
some adjustments in its shape. 

At z = the mass distribution becomes slightly less 
spherical towards the center with axis ratios dropping 
from ~ 0.55 at 200 kpc < r < r 2 oo to ~ 0.45 at 20 kpc, 
whereas the potential shape axis ratios remains constant 
at around 0.8. The triaxiality parameter remains firmly 
in the prolate regime (> 2/3) for both mass and poten- 
tial throughout the entire halo. The potential shape be- 



comes slightly more prolate in the inner regions. At ear- 
lier times, however, we detect some significant changes in 
the axis ratios, especially in the outer regions. At z = 0.5 
both q and s are larger by about 0.1-0.15 beyond 100 
kpc, and at z = 1 q is significantly larger than s, result- 
ing in a triaxial, as opposed to prolate, outer region be- 
yond ~ 200 kpc. The potential shape exhibits variations 
of comparable magnitude. Note that the host halo ac- 
cretes some fairly massive subhalos (M su b/Mh os t ~ 0.1) 
between z = 1 and z = 0.5. Dynamical friction causes 
these subhalos to spiral in to the center over a few orbits, 
and they lose most of their mass (> 99%) in this process. 
The associated redistribution of material probably con- 
tributes to the observed shape adjustments. 

The shape of the Via Lactea host halo is consistent 
with previous studies of the shapes of dark matter halos, 
whi ch have generally fo und them to be mostly prolate 
(e.g lAllgood et ail l2006. and references therein). Obser- 
vational studies of the shape of the Milky Way galaxy, 
however, have been much less conclusive. Disk warping 
(lOlling fc Merrifieldll2000l) and some models of the Sagit- 
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Fig. 2. — Host halo shape parameters for the mass distribution 
(top panel) and the potential (bottom panel) as a function of proper 
distance from the halo center at z = (solid lines), z = 0.5 (dashed 
lines), and 2 = 1 (dotted lines). Intermediate to major axis ratio q 
(red), minor to major axis ratio s (blue), and triaxiality parameter 
T (green). 



tariu s tidal stream fllbata et all l2001i iMajewski et al.l 
120031: iMartmez-Delgado et all |2004[ ) suggest that the 
Milky Way halo is close to spherical (q w s > 0.8) and 
oblate, whereas some studies of the leading arm of the 
Sagittarius stream favor a p rolate shape with s = 0.6 
rtHelmill2004t lLaw et al.ll2005h . 

The collisionless nature of our Via Lactea simula- 
tion prohibits a direct comparison between our find- 
ings and the observational constraints on the shape 
of the Milky Way host halo. Previous hydrody- 
namical numerical studies of galaxy formation have 
found that the cooling of baryons leads to significantly 
rounder halos (iKazantzidis et alj l2004t ISpringel et al.1 
120041 iBailin fc Steinmetdl2005l Tespeciallv in the central 
regions where s can increase by as much as 0.2 to 0.3. At 
the moment it is unclear how much of this effect is due 
to a potential overcooling problem in these simulations. 

4. SUBHALO SHAPES 

We now turn to the analysis of the shapes of Via 
Lactea's subhalos. Unlike the host halo, which is re- 
solved by over ten million particles even within O.l^ooi 
the subhalo shape determination is limited by numerical 
resolution. We restrict our analysis to halos containing 
at least 5,000 particles within their tidal radius, r t . As 
in paper I, we define r t as the point where the subhalo's 
spherically averaged density profile drops to twice the 
local underlying matter density. Via Lactea contains 97 
such halos within r2oo, 309 within 3r2oo, and 829 in total. 

One complication arises from the difficulty in distin- 
guishing between actual subhalo particles and those be- 



longing to the underlying host halo. Our subhalo finder 
(6DFOF, cf. paper I) assigns all particles within the 
tidal radius to the subhalo, without removing unbound 
particles. Some contribution to and Kij would thus 
come from background particles. This contribution is 
generally negligible for Zy, since the background parti- 
cles are more or less uniformly distributed throughout 
the subhalo. The host halo particles, however, signifi- 
cantly distort Kij, because their velocities are typically 
much larger and more anisotropic than the subhalo par- 
ticle velocities, especially for subhalos close to the host 
halo center. For this reason we have cleaned the sub- 
halos from background host halo particles by comparing 
the members of each subhalo to those of its progenitor in 
the penultimate simulation output at z = 0.005 (~ 68.5 
Myr before z = 0), and retained only those particles 
that also appear within r t at the earlier time. This ef- 
fectively removes all host halo particles from the subhalo 
and allows an accurate determination of K^. We have 
confirmed that the shape of Uj is not affected by this 
cleaning procedure. 

The subhalo tidal radius shrinks as subhalos pass 
through the inner regions of the host halo due to the 
increasing background density. However, as we showed 
in paper II, not all particles beyond this reduced r t are 
actually stripped from the subhalo, and r t re-expands 
as the subhalo begins its climb out of the host halo po- 
tential. Measuring the shape for all particles within the 
tidal radius thus probes more central regions for sub- 
halos closer to the halo center than for subhalos in the 
outer regions. To avoid this bias we focus on the subhalo 
shapes for all particles within the radius of the peak of 
the circular velocity curve, rvmax- This choice comes at 
the cost of a reduced number of particles, but rymax does 
not temporarily decrease at pericenter passage. 

4.1. Shape vs. Radius 

We first discuss the dependence of subhalo shape on 
the radius at which they are measured, see Fig. [3] For 
this analysis we included the complete sample of 829 ha- 
los with more than 5000 particles within their tidal ra- 
dius. Most of these halos currently lie outside of r2oo, 
but a significant fraction (78% within 2r 2 oo) have passed 
through the host halo at some earlier time (paper II) 
and have experienced tidal interactions. Similarly to the 
Via Lactea host halo, we find that the subhalos' mass 
distributions become slightly less spherical in the inner 
regions; (q) decreases from about 0.87 at 3rv m ax to 0.68 
at 0.2rv m ax, the innermost point at which we determine 
subhalo shapes. Similarly, (s) decreases from 0.72 to 
about 0.45. The 68% scatter is about 0.2 for both q and s. 
As expected the potential shape shows even less depen- 
dence on radius and is almost spherical, with (q) m 0.9 
and (s) w 0.85. The potential axis parameters also have 
a smaller 68% scatter of about 0.1. 

Unlike isolated dark matter halos, which tend to be 
prolate, we find that subhalos arc generally triaxial. (T) 
shows a slight decreasing trend with radius, but remains 
in the triaxial regime for both mass and potential shape 
throughout the range of radii that we probed. Note that 
the 68% scatter around (T) extends into both the prolate 
and oblate regime. At rvmax about 25% of all (sub)halos 
arc prolate and an equal number oblate. 
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Fig. 3. — Mean (sub)halo shape parameters vs. the major axis ra- 
dius of the ellipsoidal window used in the shape determination. We 
plot the shape of the moment of inertia tensor Iij (solid lines, filled 
symbols) and the shape of the kinetic energy tensor Kij (dashed 
lines, open symbols). The 68% scatter around the mean is indi- 
cated by the shaded region for Iij and by the thin dashed lines 
for Kij. Top panel: intermediate to major axis ratio, q = b/a. 
Middle panel: minor to major axis ratio s = c/a. Bottom panel: 
triaxiality parameter, T = (1 — g 2 )/(l — s 2 ). 

4.2. Shape vs. Distance 

Next we take a look at the dependence of sub-halo 
shapes on the distance from the halo center. For this 
purpose we use the shapes determined from all particles 
within rvmax- Fig. 3] shows that the shape of the mass 
distribution is close to independent of distance, with a 
weak but significant trend towards slightly larger axis ra- 
tios closer to the host halo center. Within r2oo the mean 
axis ratios are (q) = 0.87 and (s) = 0.74, and outside 
of r 2 oo they are (q) = 0.81 and (s) = 0.64. The poten- 
tial shape is independent of distance, with (q) = 0.93, 
(s) = 0.86. The ellipsoids remains predominantly triax- 
ial over the whole range of distance, for both mass and 
potential. 

One might have expected a stronger distance depen- 
dence, since the subhalo shapes are affected by tidal in- 
teractions, which are stronger at the center. This is likely 
a consequence of measuring the shapes within rvmax, 
which is deep inside the halo mass distribution, and 
not as strongly affected by tides. Measuring the shapes 
within the tidal radius should give a more pronounced ef- 
fect, and indeed, when we look at the trend vs. distance 
for shapes measured within r t (dashed lines and open tri- 
angles in Fig. |4|) we find that for the innermost bin the 
mean axis ratios decrease. At least part of this trend, 
however, may also be explained by the temporary de- 
crease of r t as the subhalo passes through pcricenter and 
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Fig. 4. — Mean (sub)halo shape parameters measured within 
r Vmax as a function of distance from the halo center. In addition 
to the mass shape (solid lines, filled circles) and the potential 
shape (dashed lines, open circles), in the top two panels we also 
plot the shape parameters measured within rt (dash-dotted lines, 
open triangles). The thin solid line indicates r200 = 388 kpc. 

the shape-radius dependence discussed in Section |4~T1 

4.3. Shape vs. V max 

Lastly, we consider the shape parameters as a function 
of the subhalo's Knax, a proxy for mass, Fi g. [51 Again we 
use a ll 829 (sub)halos in the simulation. (|Allgood et al.l 
2006]) found that the axis ratios of isolated field halos 
measured at 0.3r v i r decreased with increasing halo mass, 
i.e. less massive halos are more spherical. In contrast, 
we find no such trend for subhalos. With the exception 
of the highest mass bin, both q and s appear to be inde- 
pendent Of Knax- 

The highest V max bin (centered on 35kms _1 , corre- 
sponding to a tidal mass of ~ 10 10 M Q ) has (s) = 
0.57 ± 0.05, a bit les s than would be expe cted from an 
extrapolation of the lAllgood et al.l (|2006f ) (s) — Mhaio 
relation which gives (s) « 0.69 ± 0.04 at 10 10 M Q for 
(is = 0.74. All but one of the six halos in this bin lie 
within 2r2oo and ca n thus be considered subhalos. In 
lAllgood et all (|2006h shapes were measured at 0.3r V j r , 
whereas our subhalo shapes are measured at rvmax- We 
find it difficult to assign a "virial" radius to subhalos (pa- 
per II), and instead use the tidal radius r t as the outer 
"edge" of the subhalo. The median ratio of rvmax to r t 
in this bin is 0.16, and for a 10 10 M Q field halo with a 
concentration of 15 rymaxAvir = 2.163/c = 0.14, so it is 
likely that we are simply measuring the shapes farther 
in, where axis ratios tend to be smaller (see Fig. [3]). For 
comparison, (s) = 0.70 when measured at r t . 

At any rate, our mean subhalo axis ratios become indc- 
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Fig. 5. — Subhalo shape parameters as a function of the subhalos 
peak circular velocity, a proxy for mass. Panels and symbols as in 
Fig.H 

pendent of mass at lower Vmax- It is possible that tidal 
interactions cause this flattening of the subhalo shape- 
mass relationship. Another possibility is that halo shapes 
(for both subhalos and field halos) arc in fact intrinsically 
independent of mass when measured at a fixed physical 
scale, such as rvmax- In this case the mass dependence 
found at 0.3r v i r would in effect just be combination of 
radius dependence and the halo mass - concentration 
relation: 0.3r v ; r is smaller than rymax in massive, low 
concentration halos. This means that relative to rv max , 
the shapes of more massive halos are probed at smaller 
radii, where axis ratios tend to be smaller. Note that 
in this picture also t he redshift dependence observed by 
lAllgood et alJ (|2006f ) might in part be a consequence of a 
window (0.3r v i r ) which becomes larger with time due to 
its co-moving definition. Further investigations will be 
necessary to fully address this issue. 

5. SUBHALO DISTRIBUTION AND SHAPE ALIGNMENT 

Having discussed the dependence of subhalo shapes on 
their properties, we now consider the spatial distribution 
and alignment within the host halo. 

5.1. Spatial Distribution 

The spatial distribution of subhalos depends sensi- 
tively on the sample selection criterion. In paper II 
we showed that the z = distribution of Via Lactea 
subhalos is anti-biased compared to the mass distribu- 
tion of the host halo. For subhalos selected to have a 
z = mass greater than 4 x 10 6 M Q , the ratio of sub- 
halo number density to host halo mass density is sim- 
ply proportional to radius, whereas for subhalos selected 



TABLE 1 

The 3D spatial distribution of subhalos 
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Note. — The shapes of the subhalo spatial distribution for 
different subhalo samples (described in the text). is the number 
of subhalos in the sample, q and s the axis ratios, and a ,b,c the 
angle between the major, intermediate, and minor principal axis of 
the subhalo and the host halo ellipsoid. 

to have a z = V max greater than 5kms _1 , this ra- 
tio scales as the enclosed mass. Previous studies have 
found that when subhalos are selected by their mass 
or Vmax at the time of accretion, this bias is substan- 
tially reduced, and the resulting subhalo number density 
profil e more closely traces the underlying mass distribu- 
tion (|Nagai fc Kravtsov 2005; Faltenbacher fc Diemandl 
[200l . 

Here we extend this analysis to a full 3D spatial distri- 
bution of the subhalos. In addition to the two samples 
considered in paper II, we include a range of samples se- 
lected by the highest V m ax a subhalo reaches throughout 
its lifetime, a quantity we refer to as 'Vm ax , P '. Three of 
these samples are defined by lower limits, V maXjP > 5, 
10, and 15 kms -1 , and three more by including the 100, 
500, 1000 subhalos with the largest V max . p . This type 
of selection is designed to remove the bias introduced by 
tidal interactions, since the selection is performed on sub- 
halo properties prior to their interaction with the host. 
For each sample of subhalos we diagonalize a weighted 
"moment of inertia" tensor (Eq[T|) constructed from the 
z = positions of the subhalos, without regard for their 
masses. We then calculate the axis ratios of the result- 
ing ellipsoids and the angles between their principal axes 
and the principal axes of the host halo's mass distribu- 
tion measured at r2oo- 

The spatial distribution of all subhalo samples is tri- 
axial, with s comparable to the underlying host halo 
density distribution and slightly larger q. Not all of 
the ellipsoids, however, are aligned with the host halo. 
The major axis of the ellipsoid defined by all subhalos 
with M t > 4 x 1O 6 M (V max > 5kms _1 ) is tilted by 
10.1° (13.5°) with respect to the host halo's major axis. 
In general, the subhalo samples selected by Vn ax ,p are 
more closely aligned with the host halo. The best align- 
ment is found for the sample consisting of the 500 sub- 
halos with the largest V max , P , whose principal axes are 
tilted by less than 5° from the host halo's. The results 
for all sample s are s ummarized in Table [Tl Note that 
Zentner et ail (120051 ). ILTbeskind et all (|2005l . 2007), and 
Kang et all (|2007h also find that the substructure distri- 
bution is well aligned with the host halo orientation. 

The radial dependence of the subhalo number density 
is presented in Fig. [6l where we plot the ratio of the 
subhalo number density to the host halo's mass density 
as a function of r e , the ellipsoidal radius @. The ratio 
has been normalized to unity at r e = ?*200i in order to 
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r e [kpc] 

Fig. 6. — The ratio of the subhalo number density n su b(r e ) to 
the host halo mass density p(r e ) as a function of ellipsoidal radius, 
r e = ^/(x 1 ) 2 + (y'/q) 2 + (z'/s) 2 . The ratio has been normalized 
to unity at r e = T200 ■ The solid line indicates ro au b / p oc r and the 
dotted line n su b/p oc M(< r). 

highlight the radial dependence in the interior of the host 
halo. The top panel clearly shows the radial dependence 
published in paper II hold also for ellipsoidal binning: 
subhalo samples selected according to their present mass 
or Knax are anti-biased with respect to the host halo 
mass distribution, with the former scaling as n S ub/p oc r 
down to ~ 60 kpc and the latter as n su b/p oc M(< r) 
for the entire range of radii probed. The bottom panel 
shows that a selection based on V max , p removes a lot of 
this anti-bias. All six samples have a radial number den- 
sity dependence closer to the host halo mass distribution 
than the samples selected according to z = properties. 
The more restrictive the selection criterion is, the larger 
risub/ p becomes. For some of the selections the ratio even 
exceeds unity from 80 kpc to r2oo- indicating that these 
samples are spatially biased with respect to the host halo, 
i.e. their abundances falls off with radius faster than the 
underlying density distribution. This bias is probably a 
consequence of dynamical friction. The subhalos with 
the largest Knax,p were massive enough to experience 
some degree of dynamical friction and spiraled in towards 
the center. Along the way they lost mass, reducing the 
dynamical friction and preventing them from completely 
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Fig. 7. — The distribution of the cosine of the angle between the 
direction towards the halo center and the subhalos' major (solid 
lines), intermediate (dotted lines), and minor (dashed lines) axis, 
for all subhalos within r2oo- The data points indicate the fraction 
of subhalos within each bin of width 0.1. Top panel: 97 Via Lactea 
subhalos within T200- Middle panel: 520 Via Lactea subhalos out- 
side of r2oo- Bottom panel: 69 "D12" cluster subhalos within r vlT . 
The intermediate and minor axis data points have been offset by 
±0.01 in the horizontal direction for clarity. 

merging with the host. This process would preferentially 
remove such subhalos from the outer regions. 

The "top 500 V maX:P " sample comes closest to tracing 
p(r) (both in radial dependence and in the orientation 
of the shape ellipsoid), and we have plotted it in the top 
panel for direct comparison with the samples discussed 
in paper II. The decrease in n sn b/p below r e — 70 kpc is 
probably at least partially due to numerical resolution, 
since the decrease is smaller for larger, and therefore bet- 
ter resolved subhalos. 

5.2. Radial Alignment 

Via Lactea's very high numerical resolution allows us 
to investigate the orientation of dark matter subhalos 
with respect to the host halo center. We find strong 
evidence for preferential radial alignment of the subhalo 
triaxial mass distribution. 

In the top panel of Fig. [7] we plot the distribution of 
| cos6>|, the absolute value of the cosine of the angle be- 
tween each of the principal axes of the subhalo's triaxial 
ellipsoid and the radius vector from the host halo center, 
for the 97 Via Lactea subhalos within r2oo • If the subhalo 
ellipsoids were randomly distributed with respect to the 
halo center, this distribution would be flat. Instead we 
see that the major axis distribution is strongly peaked 
towards large values of | cos#|, indicating that the major 
axis preferentially points towards the halo center. Corre- 
spondingly, the intermediate and minor axes are biased 
towards low values of | cos6>|. 

In the middle panel of Fig. [7] we plot the same distri- 
bution for all halos outside of 3r2oo- The distribution 
is flat, showing no evidence for alignment of any of the 
principal axes. The alignment effect appears to only be 
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Fig. 8. — The dependence of the subhalo alignment towards 
the host halo center as function of subhalo radius at which the 
shape is measured. Only subhalos within T200 have been included. 
Solid line: the mean of |cos0| and its uncertainty. Dotted line: 
the fraction of subhalos with |cos0| > 0.8 and its Poisson errors. 
Radii as a fraction of rt (filled circles, lower axis) and as a fraction 
of 7*vmax (open squares, upper axis). 

present for subhalos physically associated with the Via 
Lactea host halo, and this is one piece of evidence for a 
tidal origin of this radial alignment. 

For comparison we have plotted (lower panel Fig. [7J 
the same distribution for the subhalos of a galaxy 
cluster scale dark matter halo, the 14 millio n part i- 
cle "D12" cluster discussed in iDiemand et all (|2004h . 
The same subhalo shape al i gnmen t is present here too. 
Ra gone-Figueroa fc Plionisl (|2007f ) and Faltenbacher et 
al. (2007, in preparation) also observe such an align- 
ment in a much larger sample of lower resolution group 
and cluster scale dark matter halos and their subhalo 
population. 

If tidal interactions are indeed responsible for the ra- 
dial alignment, then it may be expected that the align- 
ment would be more pronounced in the outer regions of 
the subhalo, where tidal effects are strongest. Indeed 
we observe just such a trend when we look at the mean 
of | cos 6* and /(|cos#| > 0.8), the fraction of subhalos 
with | cos 6*| > 0.8, versus subhalo radius, see Fig. [51 For 
this analysis we restrict ourselves to subhalos within r2oo, 
and consider radii as a fraction of both r t and rv max . 
When measured at r t , about 55% of all subhalos have 
a major axis that is aligned to within 35° of the direc- 
tion towards the halo center and (|cos0|) = 0.75. Both 
(| cos#|) and /(| cos#| > 0.8) decrease monotonically to- 
wards smaller r/r t . The mean tidal radius for the 97 sub- 
halos in this sample is (r t ) = 10.6 kpc, and (rv max ) =2.31 
kpc. Not surprisingly, the alignment signal is less pro- 
nounced when shapes are measured at rymax- However, 
even in this case the alignment is significant: (| cos#[) = 
0.58 with an uncertainty (= ^vard cos#|)/7Vs ub ) of 0.03, 
which corresponds to a ~ 2.5 ( 7 sign ificance. 

Note that while iLee et all (|2005l ) finds that the sub- 



z 

10 5 2 1 0.5 




-12 -10 -8 -6 -4 -2 
time [Gyr] 

Fig. 9. — (q) (solid line, red), (s) (dotted line, blue), and (V ma x) 
(dashed line, green) for the 19 subhalo sample (see text for a dis- 
cussion of their selection) as a function of time. The two filled 
circles at z = indicate (q) and (s) for all halos beyond 3r2oo, 
which have been less affected by tides from the host halo. The 
shapes are determined within rt- 

halo minor axes are preferentially perpendicular to the 
host halo's major axis in a study of cluster-scale subha- 
los, and Faltenbacher et al. (2007, in preparation) find a 
similar but weaker alignment signal in group-scale sub- 
halos, we find no evidence for such an alignment in the 
Via Lactea subhalo population. The minor axes of our 
subhalo population shows no correlation with any of the 
host halo's principal axes. 

6. REDSHIFT EVOLUTION 

The fact that the radial subhalo alignment is stronger 
for subhalos closer to the host halo center (see Fig. [7]) 
and that the alignment signal is more pronounced when 
the shapes are measured in the outer regions of the sub- 
halo (see Fig. [8j is suggestive of a a tidal origin of the 
alignment. In this section we present additional support 
for this hypothesis by looking at the temporal evolution 
of the shapes and orientations of a small sample of sub- 
halos, chosen according to the following criteria: 

i) the subhalos must lie within r2oo at z = 0, 

ii) they must have undergone at least three pericenter 
passages since z = 1.7, 

iii) they must have experienced significant tidal mass 
loss, AM/M = 1.0 - M t {z = 0)/M t (z = 1.7) > 0.4, 
and 

iv) they must contain more than 4000 particles at z = 
(M t {z = 0) > 8 x 10 7 M Q ). 

These constraints result in a sample of 19 well resolved 
subhalos that have experienced significant tidal interac- 
tions with the host halo. 

In Fig. Owe show the temporal evolution of (q), (s), 
and (Vmax) fo r this subsamplc. The (V max ) curve shows 
that these subhalos experience most of their mass growth 
early on and then continually lose mass due to tidal in- 
teractions with the host halo (see paper II for a more 
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TABLE 2 

SUBHALO PROPERTIES OF THE SHAPE EVOLUTION SUBSAMPLE 
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Vmax 


AM/M 


Id 
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(M ) 
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04242 


2.4 x 10 8 


15.4 


3.7 x 10 9 


28.2 


0.93 


10876 


1.6 x 10 8 


14.6 


4.6 x 10 8 


18.8 


0.66 


13351 


4.4 x 10 s 


18.8 


2.6 x 10 9 


66.9 


0.84 


13467 


1.2 x 10 s 


14.3 


2.8 x 10 9 


28.6 


0.96 


18412 


1.3 x 10 s 


12.5 


4.0 x 10 9 


32.5 


0.97 


21500 


8.4 x 10 7 


11.9 


1.5 x 10 8 


13.4 


0.45 



Fig. 10. — The proper space orbits and shapes of six subhalos (see text for a discussion of their selection) as function of time. Colors 
indicate the distance in proper kpc from the host halo center, which is indicated by the red icosahedron. The symbols depict ellipsoids of 
a constant size (major principal axis of 10 kpc) whose axis ratios and orientations are determined from all particles within the subhalo's 
tidal radius. Ellipsoids are plotted for outputs between z = 2 and z = with a stride of 4 outputs, corresponding to a time interval of 274 
Myr. The major principal axes are indicated by short solid lines. 

extensive discussion). Here we show that this tidal in- 
teraction also leads to subhalos becoming rounder with 
time. At formation they have (q) « 0.65 and (s) w 0.5, 
but these grow as they are captured and begin to feel 
the host halo's tidal field. After z = 1 (q) remains stable 
at about 0.9, but (s) is still increasing slightly, reaching 
~ 0.85 at z = 0. For comparison we also show the z — 
mean axis ratios of the 520 halos outside of 3r2oo- Ow- 
ing to their large distance from the host halo, these halos 
have experienced weaker tidal interactions. As expected 
these halos are less spherical than the tidally stripped 
sample: (q) = 0.85 and (s) = 0.72. This further supports 
the notion that tidal interactions tend to make subhalos 
rounder. 

In the following we further restrict our sample, and 
look in more detail at the orbits of individual subhalos 
and the time dependence of their shapes and alignments. 
For this purpose we hand selected five illustrative exam- 
ples of subhalos with small pericenters ranging from 9.6 
to 32 kpc. For comparison we also included one subhalo 
(Id #21500) with a higher angular momentum orbit and 
a pericenter of 69 kpc. The properties of the selected 
subhalos are summarized in Table [2j 

In Fig. [10] we present a three-dimensional visualization 
of the orbits for the six subhalos in this sample. Starting 
with the z — output and for every fourth output there- 
after (corresponding to a time interval of 274 Myr) up to 
z = 2, we have plotted at the subhalo's center of mass lo- 
cation an ellipsoid whose orientation and axis ratios are 



determined from all particles within the subhalo's tidal 
radius. The major axis of each ellipsoid is indicated by 
a short solid line. For clarity we used a fixed major axis 
length of 10 kpc, although the tidal radii vary along the 
orbit. The orbits are shown in proper coordinates in the 
rest frame of the host halo and the sides of the box range 
from -300 to +300 kpc. 

The most striking feature of these plots is that the 
radial alignment of the subhalo is preserved throughout 
the majority of its orbit. The subhalos' ellipsoids perform 
a near perfect figure rotation such that the major axis 
continually points close to the host halo's center. This 
figure rotation, however, is not seen in subhalo #21500, 
the one with the higher angular momentum orbit. Its 
orientation is almost independent of its orbital position. 

We have quantified these trends in Fig. [TTJ where we 
plot each subhalo's radius, its axis ratios q and s, and its 
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Fig. 11. — The shape parameters (measured within rt) and alignment of six subhalos (see text for a discussion of their selection) as a 
function of time. In each plot the subhalo's proper distance from the host halo is shown in the black solid line. For each subhalo the top 
panel also shows q (upper line, red) and s (lower line, green) and in the bottom panel | cos0|, the angle between the subhalo's major axis 
and the direction towards the halo center. 



| cos 9 as a function of time. The evolution of q and s 
very closely tracks the orbit of the subhalo. Every time 
a subhalo passes through pericentcr its axis ratios de- 
crease. Note that part of this decrease can be attributed 
to the fact that the tidal radii shrink during pericenter 
passage, and so we are effectively measuring the subhalo 
shapes further in, where they are intrinsically less spher- 
ical (see Section 14. ip . However, r t drops below ry max 
only for one or two outputs close to pericentcr, and for 
most of the orbit remains in the regime where the axis 
ratios arc almost independent of radius. Thus we con- 
clude that the temporary decrease in r t at pericenter is 
not sufficient to explain the full extent of the correlation 
between axis ratios and orbital position. Furthermore we 
see from the plot of | cos#| that the subhalos are indeed 
pointing close to the host halo center for the majority of 
their orbits. During the fast pericenter passages |cos0| 
drops significantly, but remains close to unity almost ev- 
erywhere else. There are a few exceptions, where | cos#| 
drops below unity even away from pericenter, for exam- 
ple at z = 0.7 for #04242 and at z = 0.4 for #10876 
and #13467, and these are most likely caused by close 
passages to other massive subhalos. 



Together these trends provide strong evidence for a 
tidal origin of the radial alignment of subhalos. As they 
orbit the host halo's center of mass, tidal forces contin- 
ually distort the subhalos' mass distribution, stretching 
them along the radial direction and compressing them 
in the perpendicular directions. This tidal distortion is 
stronger the closer the subhalo gets to the host halo cen- 
ter, but the pericenter passage occurs so quickly that the 
subhalo does not have enough time to adjust its orienta- 
tion to point towards the center. 

Our subsample, of course, was chosen to have experi- 
enced strong tidal interactions by requiring AM/M > 
0.4 and at least 3 pericenter passages. However, as we 
showed in paper II, more than half of all subhalos lose 
more than 50% of their mass from z = 1 to z = 0, and 
so tidally caused radial alignment is expected for most 
subhalos. As mentioned above, subhalo #21500 is one 
example of a subhalo that does not appear to experience 
much radial alignment. We found that of the 19 subhalos 
for which we performed a time dependent analysis, only 
three subhalos exhibit a similar lack of alignment-orbit 
correlation. 
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7. CONCLUSION 

The main conclusions of this paper can be summarized 
in the following points. 

• The shape of the Via Lactea host halo is prolate, 
and slightly less spherical (lower axis ratios) in 
the central regions. The shape of the potential 
(q w s ~ 0.8) is significantly rounder than the mass 
distribution (q w s m 0.4 — 0.55). 

• Whereas isolated halos tend to be prolate, we find 
that subhalos are predominantly triaxial. Over- 
all subhalos arc more spherical than the host halo. 
The mass shape has (q) = 0.83 and (s) = 0.68, av- 
eraged over all subhalos within 3r2oo- The poten- 
tial shape for these subhalos is very close to spher- 
ical, with (q) = 0.93 and (s) = 0.86. Like isolated 
halos, subhalos tend to be slightly less spherical in 
the central regions. 

• We find a weak trend towards larger axis ratios 
for subhalos closer to the host halo center. Within 
r 200 (?) = 0.87 and (s) = 0.74, compared with 
(q) = 0.82 and (s) = 0.65 for all subhalos outside 
of r 20 o- 

• For subhalos with V max < 30kms _1 the axis ratios 
are independent of V max . At higher V max they are 
slightly lower. 

• The spatial distribution of subhalos matches the 
prolate shape of the host halo when subhalos are 
selected by V max . p . the largest V mtlx they ever had 
during their lifetime. This type of selection also 
results in an ellipsoidal radius dependence of the 
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